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Abstract 

We consider the two-dimensional t-t' Hubbard model with the Fermi level 
being close to the van Hove singularities. The phase diagram of the model is 
discussed. In a broad energy region the self-energy at the singularity points 
has a nearly-linear energy dependence. The corresponding correction to the 
density of states is proportional to In 3 |e|. Both real- and imaginary part of 
the self-energy increase near the quantum phase transition into magnetically 
ordered or superconducting phase which implies violation of the Fermi-liquid 
behavior. The application of the results to cuprates is discussed. 
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I. INTRODUCTION 



Last decade, a possibility of a non-Fermi-liquid (NFL) behavior in two-dimensional (2D) 
systems has been a subject of many theoretical investigations. This question is especially 
interesting in connection with the high-temperature superconductors where many unconven- 
tional features, including pseudogap phenomena, are observed. Anderson has put forward 
the idea that a 2D system can demonstrate NFL behavior at arbitrary small interelectron 
repulsion U owing to a finite phase shift at the Fermi energy. Even in the absence of such 
effects, a NFL state can occur at small U because of peculiarities of electron or spin and 
charge fluctuation spectra. 

In the usual 2D Fermi liquid the imaginary part of the self-energy (electron damping) 
has the energy dependence ImS(kp,e) oc e: 2 In |e| @||, and the temperature behavior of 
resistivity is p oc T 2 In T j|| . However, these dependences do not describe experimental data 
on high-T c cooper-oxide compounds. To treat the anomalous behavior of these systems, 
Varma et al || proposed the phenomenological marginal Fermi-liquid (MFL) theory where 
Im£(kp,£) oc |e|. Then the electronic specific heat should demonstrate T InT-dependence, 
and the resistivity the T-linear behavior. A similar behavior in a broad temperature region 
can be obtained in the presence of strong antiferromagnetic fluctuations for 2D and nested 
3D systems . The linear energy dependence of the self-energy was also obtained within 
the spin-fermion theory |||| which supposes that the magnetic correlation length of the 
system is large enough. 

Another explanation of anomalous electron properties of 2D lattice systems can be found 



in the presence of the van Hove (VH) singularities [ 1Q| . In this case, the bare electron density 
of states is logarithmically divergent and to leading (second) order in U we have the behavior 
ImSikvH, s) oc \e\ In at |e| ^> (kyn is assumed to be VH point), which differs slightly 
from the linear dependence in the MFL theory. However, the important question occurs how 



this behavior changes with taking into account higher-order terms. RPA calculations [I I 



numerical calculations within the FLEX approximation |T2| and analytical analysis |12[] show 
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that not too close to van Hove filling the (quasi-) linear energy dependence of the self-energy 
takes place in a broad range of the energies. 

Within the renormalization-group approach, the behavior of the self-energy was discussed 
in Refs. [P^ , P^j . The conclusions of different approaches turned out to be contradictory: a 
possibility of a NFL state is discussed in Ref. [|13|, while a standard Fermi-liquid behavior 



was found in Ref. [14]]. It should be noted that in both the papers only part of all the 
possible interaction channels was taken into account. However, as it was discussed in Refs. 
|l5|Jl6|] , all the scattering channels are equally important, which can substantially change 
the previous results. 

In the present paper we consider the energy-dependence of the self-energy with account 
of all the scattering channels. We demonstrate that a violation of the Fermi-liquid picture 
takes place provided that the system is close to quantum phase transition into a magnetically 
ordered or superconducting state. Our approach can be considered as a generalization of 
that of Ref. JET] to the case of a non-nested Fermi surface. Note that our RG analysis takes 
into account the contribution of the vicinity of VH points only. For fillings close to VH 
one, this contribution is most singular and vertices with the corresponding momenta are 
most divergent. This can be checked by analysis of the results of complete RG approach 
which takes into account the contribution of entire Fermi surface . Unlike this approach, 
however, present approach gives the possibility to take into account all channels of electron 
scattering, as discussed in details in Ref. JT5| . Recent analysis |16[ which permits to consider 
the interplay of different channels gives the results close to that obtained earliar within two- 
patch and parquet equations approach Jl5[]. Thus we can expect that the approach based on 
consideration of the vicinities of VH singularities is reliable and gives qualitatively correct 
results for fillings close to VH one. The behavior of some physical properties and application 
of the results obtained to cuprates is discussed in the Conclusion. 
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II. THE MODEL AND QUANTUM PHASE TRANSITIONS 



We consider the t-t' Hubbard model on the square lattice: 

k i 

with the electron spectrum 

£k = — 2t(cos k x + cos k y ) — 4t'(cos k x cos k y + 1) — [A (2) 

Hereafter we assume t > 0, t' < (which is the case for hole-doped systems), < \t'/t\ < 1/2. 
For t' = the Fermi surface is nested, which results in peculiarities of physical properties 
17,^1 ■ However, nesting is removed for t'/t ^ 0. For arbitrary t' /t, the spectrum @ contains 



VH singularities connected with the points A — (jr, 0), B — (0,7r). The chemical potential 
fi is determined by the electron concentration n and can be obtained from the condition 



n 

k 



£/k (3) 



with /k = /(^k) being the Fermi function, so that fi = corresponds to Van Hove filling. 
Being expanded near the VH singularity points, the spectrum (|2|) takes the form 

£k = -2t(sin 2 (pkl - cos 2 tpky) - n (4a) 
= 2t(cos 2 (fk 2 x - sin 2 (pk y ) - /i (4b) 

where k x = tt — k x , k y = it — k y , <p is the half of the angle between asymptotes at VH 
singularity, 2ip = cos _1 (— 2t'/t). 

We have a set of quantum phase transitions (QPT's) for the model (0), see Ref. [|T^ 



We restrict ourselves to the case of small enough U < Qt, so that we may neglect the 
correlation effects connected with Hubbard's subband formation (e.g., the Mott-Hubbard 
metal-insulator transition). Provided that the electron concentration is not too close to its 
van Hove value, i.e. the chemical potential satisfies \fi\ > fi c , fi c being the critical value, we 
have the normal (paramagnetic and non-superconducting) phase. When approaching the 



VH singularity, a QPT to antiferromagnetic, superconducting or ferromagnetic state occurs. 
The type of the ground state depends on t'/t and U. In particular, for U — At we have the 



antiferromagnetic ground state at \t'/t\ < 0.30 and ferromagnetic one for larger \t'/t\ |T5 
The critical electron concentrations for stability of these phases differ by several percents 
from van Hove filling. 

III. ELECTRON SELF-ENERGY IN THE SECOND ORDER 



Consider first zero-temperature perturbative results. To second order the expression for 
the electron self-energy has the form 

E^(k, e) = U 2 J2 /p+q " k(1 - /p - /q) + /p/q (5) 

This contribution can be represented as a sum of three diagrams (Fig. 2 a-c): 

£ (2) (k, e) = Si(k, e) + E 2 (k, e) + E 3 (k, e) (6) 

When picking out the singularities we can put e p = e p or e p = e p for p being close to A(B) 
van Hove points. Note that the term Si was investigated earlier [|13|,[14],^2|] . We restrict our 
consideration to the VH points at the Fermi surface, k^ = (0, n) or (vr,0). The calculations 
at |e| yield: 

Re£i(k F ,e) = -(ln2)(^ /sin2^) 2 eln 2 (A 2 |t/e|) 

v4 23 ln(A 2 |2t/£|) |£|<2tcos2^ 
ReZ 2>3 (k F ,e) = -(g /sin2tp) 2 e{ * (7) 

I fc 2 , 3 ln21n 2 (A 2 |2t/£|) |£|>2tcos2y? 

where go = U/(An 2 t) is the dimensionless coupling constant, A ~ 1 is the ultraviolet mo- 
mentum cutoff, k 2 = 1, k 3 = 2/3, 

cos 2<p 1/ cos 2ip 

f dx , e x — x 1 f dx , e, + x 
A 2 = I — In + - / — In 



cos 2ip 



1 f dx e x 



+1 / _i n _^_ ( 8a ) 
2 J x e x + x 



-oo 



cos 2ip 1 / cos 2ip 

, dx , / x \ 1 f dx 2(e x + x) 
A 3 = I - —In (-—)+- / — In v 



2e x + x \ 2e x J 2 J 2e x + x x 

cos 2ip 



1 f dx , x , , . 

+ - / h-7 r 8b 

— oo 

with e x — (x — cos2</?)(l — x cos 2y?)/ sin 2 2yj. For small energies |e| <C \fi\ the real part 
remains linear in energy with the difference that logarithmical divergences are cut at \fi\ 
rather than at |e|. The corresponding imaginary parts at |e| ^> read 

ImEi(k F ,e) = -(7rln2)(Vsin2^) 2 |e| ln(A 2 |t/e|) 

B 23 \e\ < 2tcos2y2 

Im£ 2 , 3 (k F ,£) = -n(g / sm2<p) 2 \s\ I "' (9) 

I fc 2 , 3 ln21n(A 2 |t/e|) |e|>2tcos2y? 

where 

B 2 = (A 21 + A 23 )6(e) + A 22 6{-e) 
B 3 = (A 31 + A 33 )9(-s) + A 32 9{e) 

and A 2n and A 3n are n-th summands in the definition of A 2 and A 3 , Eq. @, respectively, 
0(e) is the step function. At t'/t ^ we have ImE(k j p,e) 7^ ImS(ki?, — e) because of the 
absence of the particle-hole symmetry for the less singular terms X 2j3 (kF, e). This fact 
results in an asymmetry of the electron density of states near the Fermi level and can be 
important for some physical properties, e.g., thermoelectric power. Unlike the real part of 
the self-energy, the imaginary part changes its dependence to quadratic one at |e| <C 
demonstrating a conventional Fermi-liquid behavior in this region. 

IV. RENORMALIZATION 

Both real and imaginary parts of the self-energy contain large logarithms at <C |e| <C t. 
Therefore we can introduce the logarithmic variable A = lr^Ali/el 1 / 2 ). Then the leading 
terms of the expansion in the powers of interaction strength can be written down as 



S(A) = AXg 2 e(l + C l9o X + D l9 2 X + ...) 

To perform the summation of leading logarithms in the self-energy, we introduce the vertices 
7i(A), i = 1, ...4 (Fig. 2), and consider the renormalization of E. As discussed in Refs. 
|, |2^JI5| , 7i(A) can be determined from the renormalizat ion-group (RG) equations 

71 = 2g?i(A)7i(7 2 - 71) + 2d 2 7i74 - 2 <i 3 7i72 

72 = di(\)(il + 7 3 2 ) + 24(7i - 72)74 - 4(7i + 7 2 2 ) 
-2rfo(A) 73 74 + 2rf!(A)73(272 - 71) 



7a 



V 4 = -ci (A)(7 3 2 + 74) + rf 2 (7i 2 + 27i72 - 2 7a 2 + 7 4 2 ) 



(10) 



where 7- = dji/dX, 



d (X) = 2c A; 



4(A) 



A, A < 2z Q 



z q , X > 2z Q 
d 2 = 2z ; d 3 = 2c Q 



(11) 



The quantities 



z = c = l/sin(2y>) = 1/V1 - B? 



are the prelogarithmic factors in small-momentum particle-hole and particle-particle bubble, 
while 



Zq = ln[(l + VT^)/R] 



c Q = tan-^fi/Vl -R 2 )/R 



(12) 



are the prelogarithmic factors in particle-hole and particle-particle bubble with momenta 
close to Q = (71,7c), R = —2t'/t. Equations (10) should be solved with the initial condition 
7i(0) = 9o- Because of the presence of double-logarithmic terms, the corresponding RG 
equations are only approximate. However, the comparison of the results of their solution 



15 1 with the parquet approach |15| and the RG approach which takes into account the con- 



tribution of the whole Fermi-surface [IS] shows that they reproduce well the renormalization 
of the couplings. 

The magnetic or superconducting instabilities manifest in the divergence of the vertices 
7; (A) at some critical scale A c . This is connected with the critical energy scales discussed in 
Sect. II as \i c ~ T c ~ Aexp(— 2A C ), for a detailed discussion see Ref. fl5|] . (In the absence 



of interlayer coupling the quantity T c has the meaning of a temperature of crossover into 
the state with pronounced short-range order (or pseudogap state) rather than of a phase 
transition temperature.) For A close to A c the solution of the Eqs. (|10D can be represented 
in the form 

= A^A (13) 
The corresponding correlation length is given by 

r 2 = QM(A c -A)/t (14) 

where is the universal number. For T \fi\ we should stop the scaling at 

A* = lntA^Vmax^l,^) 1 / 2 ] (15) 

while T 3> at 

A* = ln[At 1/2 / max(T, |e|) 1/2 ] (16) 

The condition > /i c or T > T c guarantees that A* < A c , i.e. the system is not ordered 
and RG approach is applicable. 

Now we consider the renormalization of the self-energy. We follow the method of Refs. 
25^p7[ . Defining the scale-dependent quasiparticle residue 



z{\) = z 1 {x)z 2 {\)z 3 {\) 

where Zj (i — 1, 2, 3) is the contribution of z-th diagram, we have the RG equations 



dlnZi(A) 

dX 
d\nZ 2 (X) 

dX 
dlnZ 3 (X) 

dX 



-(81n2)A7 2 /sin 2 2^ 

-D 2 (A)( 7l 2 + 7 2 - 7l72 )/ sin 2 2<p 

-D 3 (A)73 2 /sin 2 2^ 



(17) 



Here 



£>2, 3 (A) 



4fc 2 , 3 Aln2 A < (1/2) ln(l/cos2<^) 
A 2i3 A > (1/2) ln(l/cos2y?) 

and the summation over spin indices in vertices is performed. Then the real part of the 
self-energy can be found as 



ReS(k F ,e) = e\nZ(X*) 



(18) 



(see, e.g., Ref. |27j). After calculating ReS, the imaginary part of the self-energy can be 
obtained from the Kramers-Kronig relations. 



V. RESULTS OF CALCULATIONS 



First, we demonstrate our approach in a simple case with the only non-zero vertex, 



74 7^ 0. As discussed in Ref. Hl5| , this case corresponds to t' — > —t/2, the ground state at 
van Hove filling being ferromagnetic (flat-band ferromagnetism) . Then we have 



74 



9o 



1 + g (c X 2 - 2z X) 
The vertex (|I~9"D diverges in the critical point 



(19) 



1-y/l- l/(z o9o ) 



(20) 



2 max(/i, T) 

where we have put zq = cq in A c according to (|12"D . From (|i~7|) we obtain Z 2 = Z% = 1 and 

4 In 2 £ 2 A(A-1 



Zi = exp 



sin 2 2cp (1 - g z ) (1 + 9 o^oA 2 - 2g^ A) 



(21) 



For A being close to A c we have 



7 4~£ 2 (T = 0)~(A C -A)- 1 (22) 

As well as in the nesting case (see, e.g., Ref. |T7j|), the quasiparticle weight Z vanishes 
exponentially in the QPT point. Note that this vanishing is much faster than the inverse- 
logarithmic dependence 

z ~ • " K °- 35 (23) 

obtained in Ref. ]13| where only one scattering channel was taken into account. Besides 
nearly-linear dependences in e (with logarithmic corrections), both real and imaginary part 
of the self-energy contain large e-dependent factors of the order of 1/(A C — A) 2 oc £ 4 (0), which 
occur because of the divergence of the vertex (|19D at A — > A c . 

As discussed in Ref. [|HJ, for \t'/t\ being not very close to 1/2 we have an interplay of 
all the scattering channels, so that we have to solve Eqs. (|10D, (|17|) numerically. We also 
calculate the quasiparticle spectral weight in VH points of the Fermi surface, 

ACk f] = _I hngky^g) 

1 VHl } 7r[£ + //-Re£(k y tf, £ )] 2 + [Im£(k yH)£ )] 2 1 } 

The results of the calculations for t' = — 0A5t, \ijt — 0.3, T = (nearly ferromagnetic 
ground state, fi c /t = 0.04) are shown in Fig. 3. As well as in the second order of the 
perturbation theory, at \e\ <ti we have the 2D behavior of self-energy 

ReS(k F , e) oc e, ImS(k F , e) oc e 2 , (25) 

while \fi\ <C \e\ t it changes to van Hove behavior, 

ReS(k F , e) oc e ln 2 (At/e), ImS(k F , e) oc |e| In(At/e) (26) 

Finally, at e = Eq ~ t the real part of self-energy has a maximum and then decreases with 
farther increasing of e, while imaginary part is almost a constant at £ ~ t. The slope of RG 
results in the region \\x\ <C |e| t is substantially higher that the results of the second-order 
perturbation theory, the imaginary part of the self-energy at \e\ > |/i| also becomes large 
enough. The peak in the quasiparticle weight obtained within RG approach is more broad 
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than in the second-order perturbation theory, which is the consequence of larger electron 
damping in the letter case. With the increasing the chemical potential (so that the system is 
moved away from QPT), the higher-order renormalizations become less important and the 
behavior of the self-energy reproduces the result of the second-order perturbation theory. 

The imaginary part of self-energy at finite temperature and e = (the inverse quasipar- 
ticle lifetime at the Fermi surface) can be obtained from the scaling arguments, 

7(T) = -[ImS(k F ,£)| T=0 ]^T (27) 

Thus it is given by the same Figs. 3b with the replacement e — > T and also demonstrates 
the linear behavior in some temperature region. 

In the above consideration we neglected completely the damping of the particle-particle 
and particle-hole excitations as well as other non-singular contributions. These non-singular 
terms can be neglected provided that the condition 

£~ 2 »max(H,|e|)/t (28) 

is satisfied, i.e. not too close to quantum phase transition into the corresponding ordered 
state and at not too high energies. Close to the quantum phase transition into the ordered 



state, the peaks in the self-energy should become asymmetric, as discussed in Ref. [11], and 



the pseudogap can be formed. On the other hand, in the limit of large enough i.e. at 

max(|e|/t,£- 2 ) < \lA/t (29) 
the spin-fermion theory [0] is applicable. 



VI. DISCUSSION AND CONCLUSIONS 

In the present paper we have investigated within a scaling approach the energy depen- 
dence of the real and imaginary part of the electron self-energy in the presence of van Hove 
singularities. We have restricted ourselves to the regions of the fi — T phase diagram, where 
the Fermi level is not too close to VH points (disordered ground state, \i c < \fi\ <C t with 
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fi c is the critical chemical potential) or temperature is above the critical value (T > T c with 
T c ~ fie). 

Provided the system is not too close to QPT (at fi c < \fi\ <C t), the imaginary part 
of the self-energy at VH points demonstrates in a broad energy region <C |e| C i a 
nearly-linear behavior, ImE(kjr,e) oc |e| ln(t/|e|), which is close to that in the marginal 
Fermi-liquid concept || (however, in MFL the linear behavior of ImE(ki?,e) takes place for 
arbitrary kp). The real part of the self-energy behaves as eln 2 (t/|e|). The linear energy 
dependence of the self-energy (^) and the behaviour of the slope with changing \i are in 
agreement with the RG results of Ref. |19| which takes into account of the contribution of 
the whole Fermi-surface. These dependences are also similar to those obtained within the 
spin-fermion model for a nearly antiferromagnetic state |8]|§, although in our case they have 
a different nature and are governed by van Hove singularities themselves rather than by 
closeness to antiferromagnetic quantum phase transition. The role of a characteristic spin- 
fluctuation frequency ui s f, which separates the Fermi-liquid and MFL regimes, belongs in 
our case to the chemical potential (// = corresponds to VH filling). Another difference 
is that in the presence of VH singularities the linear dependence of self-energy takes place 
already in the weak-coupling regime. 

Near QPT, the renormalizations become important because of the large ground-state 
correlation length which enters renormalized vertices. Therefore both real and imaginary 
parts of the self-energy increase considerably as £ 2 (T = 0). Such an anomalous behavior also 
implies a strong violation of the Fermi-liquid and even MFL picture. Note that the anomalies 
under consideration may induce the electron topological transition with the truncation of 
the Fermi surface (see, e.g. Ref. ||28|| ). The RG approach used is not able to describe the 
magnetically ordered or superconducting state. By this reason, the renormalized-classical 



regime T < T c , < /i c should be considered within other approaches, see, e.g., Ref. |29|j . 

In 2D situation, the physical properities near QPT should demonstrate singularities 
which are stronger than those in the 3D case. The correction to the electron density of 
states N(e) reads 
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SN & = - E 



kcr 



ReS(k,e)5'(e - e k ) + -Im£(k,e)/(e - £ k ) 2 

7T 



(30) 



Main contribution to the integral comes from the vicinity of VH points where the bare 
density of states is logarithmically divergent. Taking into account that ReE^k, e) oc 
£ln 2 max{|£|, |e k |}, we obtain for the first (coherent) term in the square brackets which 
originates from renormalization of quasiparticle spectrum 

5N coh (e) (x\n 3 (t/\e\),e^> fi. (31) 

The calculation of the second (incoherent, non-quasiparticle) term requires the full form of 
E( 2 )(k,e), Eq.(||), and leads to the result 

8N incoh (e) ocln 2 (t/|£ + /i|). (32) 

Although this divergence is slightly weaker than of the coherent term, it is not cut at 
e = — \i. Thereby the bare VH singularity becomes considerably enhanced. Note that the 
divergence of the density of states together with its asymmetry in e may lead to peculiarities 
of thermoelectric power owing to impurity scattering, cf. Ref. [J7]. 

To leading (second) order the contribution to electronic specific heat owing to VH sin- 
gularities has the form 8C oc Tln 3 (t/ max{|/i|, T}). The resistivity (inverse transport relax- 
ation time) should demonstrate at T > /j, the behavior p oc Tln 2 (t/T). The calculations are 
similar to those of Ref. @ for the antiferromagnetic state, extra logarithmic factors coming 
from VH singularities. The crossover from quadratic to nearly linear temperature depen- 
dence of resistivity is confirmed by experimental data for cuprates (see, e.g., the results of 



Ref. [§0J for the LaSrCuO system). 

Thus the divergences in the many-electron system with VH singularities are stronger 
than those in the MFL theory. Near QPT, we can expect that all the physical properties 
are strongly renormalized, the renormalizations being dependent on the type of the ordered 
phase. This problem will be considered elsewhere. 

Finally we consider the application of the results obtained to cuprate systems. A nearly 
linear energy dependence of Im£(k>, e) at VH points, which is similar to our results, was 
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observed for the system Bi2212 in ARPES experiments |3T]. For the system La2_ x Sr a; Cu04 
with the doping x c ~ 0.2 (which is slightly larger than the optimal one) the Fermi surface 
crosses VH points J32|. The density of states specific heat coefficient and Pauli sus- 
ceptibility ]53] substantially grow near this doping. The mass enhancement factor m*/m 
demonstrates a similar behavior. The additional experimental investigations of the self- 
energy near (ir, 0) point of La2_ a; Sr a ;Cu04 would be interesting in this respect. 
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Figures 

T 
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Fig.l. A qualitative T — \l phase diagram in the vicinity of quantum phase transitions. 
The chemical potential \i is referred to the Van Hove singularity. Bold line denotes the 
ordered ground state. 




Fig. 2. The second-order diagrams for the electron self-energy. Solid and dashed lines 
correspond to electrons with momenta close to (0, n) and (tt, 0) van Hove singularities re- 
spectively 
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Fig. 3. The energy dependences of real (a) and imaginary (b) parts of the self-energy at 
VH points, and the spectral weight (c) for t'/t = -0.45, U = At, n/(2t) = 0.15 and T = 0. 
The dashed line corresponds to the second-order perturbation result. 
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